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Abstract 

In the frame of the scalar field model with non minimal kinetic coupling to 
gravity, we study the cosmological solutions of the Chaplygin gas model of 
dark energy. By appropriately restricting the potential, we found the scalar 
field, the potential and coupling giving rise to the Chaplygin gas solution. Ex- 
tensions to the generalized and modified Chaplygin gas have been made. 



1 Introduction 

A wide range of cosmological observations indicate that the universe has entered a 
phase of accelerating expansion, which becomes one of the important puzzles of the 
contemporary physics. Those observations include the type la supernovae (Snla) 
standard candles pQ, [2J, the angular location of the first peak in the CMB power 
spectrum [3] and baryon acoustic oscillations of the matter density power spectrum 
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[I] . These evidences represent a great stimulus for theoretical work and originated the 
concept of Dark Energy. The Dark Energy (DE) Models attribute the observed ac- 
celerating expansion to an unknown energy component with negative pressure, which 
dominates the universe at recent cosmological times. The data can be accommodated 
with a high degree of accuracy, in the ACDM model which combines the baryons with 
conventional cold dark matter (CDM) candidates and the cosmological constant A, 
which accounts for the dark energy. However a non zero cosmological constant raises 
the coincidence problem (why the DM and DE are comparable today), and is plagued 
by fine tuning problems [5], [6], [7]. Alternatively, dark energy may be described by 
different dynamical scalar field models with time-dependent equation of state, like 
quintessence [8], [9]; string theory fundamental scalar known as tachyon [10J; K- 
essence models involving a generalized form of the kinetic energy |llj.|12]: scalar field 
with negative kinetic term, which provides a solution known as phantom dark energy 
[T5] (see [H] for a review). An alternative description of DE may be given by perfect 
fluids with adequate equation of state, like Chaplygin gas [15], [16]. The above men- 
tioned scalar field models assume that the DM and DE are of different nature, while 
in the Chaplygin gas the DM and DE appear as manifestations of the same single 
fluid at different cosmological epochs, showing a dust-like matter behavior at early 
times and behaving as cosmological constant at late times [15], [16J. 
Another description of DE is provided by the scalar-tensor theories which contain a 
direct coupling of the scalar field to the curvature, with the advantage of giving a 
mechanism to evade the coincidence problem, and naturally allowing (in some cases) 
the crossing of the phantom barrier [17] . These theories with different couplings to 
the curvature appear as low energy limit of several higher dimensional theories, and 
provide a possible approach to quantum gravity from a perturbative point of view 
|18j . A coupling between curvature and kinetic terms appears as part of the Weyl 
anomaly in = 4 conformal supergravity [IHIEO]. A model with non- minimal deriva- 
tive couplings was proposed in [21] , [22] , [23] in the context of inflationary cosmology, 
and recently, non-minimal derivative coupling of the Higgs field was considered in 
|24j . also as inflationary model. In |25J a derivative coupling to Ricci tensor has been 
considered to study cosmological restrictions on the coupling parameter, and the role 
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of this coupling during inflation. Some asymptotical solutions for a non-minimal ki- 
netic coupling to scalar and Ricci curvatures were found in [2S], and quintessence 
and phantom cosmological scenarios with non-minimal derivative coupling have been 
studied in [27]. A scalar field with kinetic term coupled to a product of Einstein 
tensors has been considered in [28J. Non-minimal coupling of scalar fields (including 
kinetic terms) with modified f(R) theories have been also considered to solve the DE 
problem in [29], [30], [31], [32], [33]. 

In this paper we consider an explicit coupling between the scalar field the kinetic term 
and the curvature [MJ [351 ES], as the source of DE and will establish the connection 
with the Chaplygin gas, by obtaining the solution of the field equations that repro- 
duce the cosmological evolution as given by the perfect fluid obeying the equation 
of state of the Chaplygin gas [16]. The Chaplygin gas has attracted much attention 
in cosmology, as it allows to interpolate between a dust dominated phase of the evo- 
lution of the Universe in the past, and an accelerated one at recent time. On the 
other hand, as in the case of the Chaplygin gas, the scalar field with the non-minimal 
kinetic coupling to the curvature, has been shown to be successful in the description 
of the DM and DE without introducing separately the DM term (see [351 ES]), i-e. 
it describes the DM and DE as manifestations of a common scalar field at different 
epochs. For this reason, it would be interesting to analyze the connection of the 
scalar field with non-minimal kinetic couplings with Chaplygin gas, which will be the 
subject of study in the present paper. A potential and coupling function that may 
give a dynamical description to the Chaplygin gas, have been found. The generalized 
and modified version of the Chaplygin gas have been also considered. 

2 The Models and Field Equations 

The scalar field with kinetic couplings to curvature is given by [35] : 
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And we will use the flat FRW background metric given by the line element 

ds 2 = -dt 2 + a 2 (t) [dr 2 + r 2 (d6 2 + sin 2 6d<j) 2 )] (2.2) 
where a is the scale parameter. The dimensionality of the coupling constants £ and 



1] depends on the type of function F ((/)). Taking the variation of action (2.1) with 
respect to the metric, we obtain a general expression of the form 

1 
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k 2 T, 



(2.3) 
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where k 2 = 8ttG and the tensor T^ u represents the variation of the terms which 
depend on the scalar field and can be written as 

T, u = Tl + T% + Tl (2.4) 
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where Tff u , T^, T^ v correspond to the variations of the minimally coupled terms, the 
£ and the r\ couplings respectively. Due to the interaction between the scalar field 
and the curvature, the derived expressions for the density and pressure for the scalar 
field can be regarded as effective ones. From now on, and in order to simplify the 
equations (the field equations will contain only second order derivatives) we will use 
the restriction t] = — 2£, which is equivalent to a coupling of the kinetic term to the 
Einstein tensor (see [22], [23]). Evaluating the 00 and 11 components of the Eq. 



(2.3) in the spatially-flat Friedmann-Robertson- Walker (FRW) background (2.2), it is 
obtained (with the Hubble parameter H = a/a, and for homogeneous time-depending 
scalar field) 
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where "dot" represents the derivative with respect to the cosmological time t. Taking 



variation in (2.1) with respect to the scalar field in the FRW background, gives the 
equation of motion as follows 

4> + 3H<j) + %- + 3£H 2 ( 2F(0)0 + j + l8£H 3 F((j))<j) + 12£MF(0)0 = (2.7) 
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where the first three terms describe the minimally coupled field. In what follows we 



will study cosmological solutions to Eqs. (2.6) and (2.7) giving rise to accelerated 



expansion, and according to the cosmological scenario described by the Chapkygin 
gas solutions. 

The Chaplygin gas is described by the following equation of state (EoS) 

A 



P 



(2.8) 



where p and p are respectively the pressure and density, and B is a positive constant. 
The continuity equation takes the form 

A" 

P. 

This equation can be easily integrated in the variable x = log a, yielding 
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which according to the Friedmann equation gives the following Hubble parameter 
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This solution has the known advantages of describing the presureless matter domi- 
nance stage (epoch) of the universe at early times (a << 1, normalizing the current 
value of a to 1), and the future universe dominated by the cosmological constant, 
entering in a de Sitter phase at a » 1. In the next section we will consider this 



solution to integrate the equations (2.5) and (2.7) with respect to the scalar field 



and the coupling F, and in this manner we will obtain a description of the Chaplygin 
gas cosmology in the frame of the scalar field with non-minimal kinetic coupling to 
curvature. Then, we will follow the same procedure with the generalized and modified 
Chaplygin gas models. 



3 Standard, generalized and modified Chaplygin 
gas solutions 



In order to integrate the Eqs. (2.5) and (2.7) for a given Hubble parameter, we should 



impose additional restrictions on the scalar field potential in order to consistently find 
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the rest of the variables, as follows (see [35]). 

In terms of the variable x = log a, and defining the function 6(x) = <f)' 2 , the Eq. (2.7) 
can be written as (after multiplying by <p) 

I A. (h 2 9) + 3H 2 9 + ^ + %H 2< ^F9 + 3£# 4 -^-(F#) + l&^FO = (3.1) 



From Eq. (2.5), changing to the variable x, we can write the product F<ft' 2 = F9 as 
following 
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taking the derivative of Eq. (3.2) and replacing F9 and d{F9)/dx into Eq. (3.1), we 
arrive at the following equation involving 9, H and V 
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In this manner, we obtain a first order differential equation for the functions 9, H and 



V. In order to integrate the equation (3.3), and thanks to the fact that the functions 
9 and V are separated, we can impose a restriction on the scalar field potential given 
by the equation 
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which simplifies the Eq. (3.3): 
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The Chaplygin gas solution 



In order to consistently solve the Eqs. (3.4) and (3.5), we propose the expression 
for the Hubble parameter H 2 , given by that of the Chaplygin gas (2.11 ) 
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Defining the scaled Hubble parameter H = H/Hq, Eq. (3.6) can be written as 

1/2 



E 1 



A + Be 
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where the parameters A and B are now dimensionless and are given by 
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where A and B satisfy the flatness condition (considering the Chaplygin gas domi- 



nance 



A + B = 1 



(3.9) 



Replacing H 2 in (3.4), and defining the dimensionless scalar potential V = k 2 V/Hq, 



after changing the Eq. (3.4) to the "tilde" variables and integration we obtain the 
scalar field potential 

1 5 Ae 6x 
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where C is the integration constant. Replacing H 2 in (3.5) we get the following 
expression for 9 

e(x) = <p' 2 = 6oe ~ 6 \ 1/4 , (3.ii) 

(A + Be- 6x ) 

where 8q is the integration constant. Integrating the square root of this last equation, 
we obtain the scalar field as (considering the (— ) sign root) 
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Finally, the coupling function F is found by replacing the Eqs. (3.7)-(3.11) in the 



Friedmann Eq. (|3.2|), giving the result 

k 2 6< 



Fx) 



3^ 2 H 2 6 



g[x 



1/4 ^ t/Q 6x C „3x 

6 3 



e ~+ e 6*o (a .)i/4 2Fl 



5 



1 5 



6,r 



'2'4' 5 

(3.13) 

where g(x) = A + Be~ Gx . Although we can not have an analytic expression for the 
potential in terms of the scalar field, we can illustrate the behavior of the potential as 



7 



showed in fig.l, for a given set of parameters. The constant C is selected so that the 
time variation of the gravitational coupling does not exceed the observational limits 
(see below). As can be seen from Fig. 1, for the selected values of the parameters, 
the potential is a monotonic decreasing function of the scalar field. The decreasing 
runaway behavior of the potential describing dark energy, is a key fact for a realistic 
cosmological model 
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Fig. 1 The evolution of the potential for the Chaplygin gas solution with the scalar 

field for A = 0.7, C = 11. 

The generalized Chaplygin gas solution 

A generalization of the Chaplygin gas has also been considered to describe the dark 
matter and dark energy. It's equation of state is given by 



Pg = -- 



Pg 



(3.14) 



where pc and pc are the energy density and pressure of the generalized Chaplygin 
gas, A is a positive constant and a is considered to lie in the range < a < 1, which 
guarantees the stability and causality (see below) [T5], [IB], [35]. Note that a = 1 



corresponds to the original Chaplygin gas. The EoS (3.14) has an equivalent field 



theory representation in a generalization of the Born-Infeld theory [39] (in the scalar 
field representation used in [39], the Born-Infeld Lagrangian density is reproduced for 
a = 1). solving the continuity equation (p G + 3H(p G +Pg) = 0), leads to the energy 
density in terms of x — log a 



p G = (A + Be- 3{a+1)x ) 1+0 



(3.15) 



where B is a positive integration constant. This density for a > — 1 describes the 
matter dominated phase at a << 1 and at the limit a » 1 describes the de Sitter 



phase dominated by constant density. At the limit a — > 0, Eq. (3.15) reproduces 
the ACDM model. This model also captures the attention because of it's connection 



with string theory and supersymmetry: the EoS (3.14) can be obtained in the Nambu- 

and has 



Goto action for ci-branes moving in a (d + 2) dimensional space time 
supersymmetric generalization [41J. As a criteria to constraint the constant a we can 



use the sound speed for the fluid described by Eq. (3.14), given by 

2 dpc ctA aA 

Cs = d^ = A + Be- 3 ( 1+a > = A + B{1 + z) 3 U+«) 

where the last is written in terms of the redshift z (e~ x = (1 + z)). At future z — > — 1, 



(3.16) 



and as follows from (3.16) c s — > a. Therefore, stability requires a > 0, and causality 
requires a < 1 (i.e. respects the speed of light limit), limiting a to the interval 
< a < 1 (note that there is not such restriction on a at high redshift, during the 
matter dominated epoch). However some authors have considered the "forbidden" 
region a > 1 if the sound speed c s is treated as group velocity [42], (see also |43] for 
superluminal sound speed). According to some observational studies, for the case of 
the pure GCG the a < 10 -4 values are favored, which is very close to the ACDM limit 
|44j . This tight restriction may be increased if we consider additionally the barionic 
matter component (a < 10~ 2 ) [45], or adding cold dark matter and barionic matter 
components (a < 0.2) [46] . On the other hand, based on observations of the barionic 
power spectrum, the GCG plus barion matter is favored for a > 3 [12], [43] . 



Replacing the density (3.15) in the Friedman equation gives the following scaled 
Hubble parameter 
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A + Be 
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where the dimensionless A and B are given by 
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(3.17) 



(3.18) 



and satisfy the flatness condition (for pure generalized Chaplygin gas content) A+B 



1. Replacing (3.17) in (3.4) and solving the Eq. (3.4) in "tilde" variables we obtain 
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the following solution for the scalar potential 
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where C is the integration constant. Solving (3.5) with H 2 given by (3.17), we find 
the expression for 9 



9{x) = <p 
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where #0 is the integration constant. Integrating the square root of (|3.20|) it follows 
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And the corresponding coupling function, as follows from (3.2) and (3.17 3.20) is 
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(3.22) 
i 

where g a (x) = i^A + Be~ 3 ^ 1+a ^ x j 1+a . In all equations the dependence on the redshift 
z or in the scale factor a is obtained by replacing e~ x = 1 + z = a -1 . In fig. 2 we 
plot the behavior of the potential in terms of the scalar field for two values of a 
corresponding to the "physical" region < a < 1, and the superluminal region a ~ 3 
which also accommodates in some astrophysical observations [42], [43] . 
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Fig. 2 The potential for the GCG versus the scalar field for 
A = 0.5, a = 0.6, C = 25.7 (dashed), and A = 0.5, a = 3, C = 3.1. iVote toe 
decreasing behavior, which is an important characteristic for dark energy potentials 

The scalar field potentials for the non-minimally coupled scalar field, that reproduce 
the dynamics of the Chaplygin and generalized Chaplygin gas are decreasing functions 
of the scalar field. Note the runaway behavior for all curves, which are characteristic 
of well behaved dark energy potentials [37] , [38] . The parameters are chosen in such 
a way that V(<f>) is definite positive, at least for z > —1. For other choices of the 
parameters, it can be shown numerically that the curves can have one maximum. 
C is used to accomplish the observational restrictions on the time variation of the 
gravitational coupling. 



So, we have reconstructed the scalar model with kinetic couplings to curvature (2.1), 
for a given Hubble parameter describing the CG and GCG cosmologies. In both cases 
we exploited the additional degree of freedom represented in the coupling function, 
to constraint the scalar potential in a way consistent with the Friedmann equations. 
A more general formulation of cosmological reconstruction method (in time and x 
variables) for a number of modified gravities including scalar tensor theories, have 
been performed in jUJ, [48] . In the scalar tensor theories considered in [UJ, |48j . the 
scalar potential and couplings are reconstructed by using an apropriate redefinition 
of the scalar field, and giving the particular type of cosmological evolution encoded 
in H. The reconstruction was considered in the cosmological time and the e-folding 
variable x, and concrete examples of accelerated late time cosmologies have been 
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provided. In this works for the case of f(R) gravity, the reconstruction was achieved 
by introducing an auxiliary scalar field. 
The modified Chaplygin gas solution 

The modified Chaplygin gas (MCG) is defined for the equation of state 

P = Bp - 4 (3-23) 

A generalized version of an equation of state that includes CG, GCG and MCG 
have been considered in and [5U]. Integrating the continuity equation gives the 
following energy density 

where = 3(1 + B)(l + a) and C is the integration constant. This equation has 
the same functional dependence on x or a, as the case of the GCG, except that ad- 
ditional limits can be obtained due to additional factor (1 + B) in the exponent of 



Eq. (3.24), with respect to the corresponding one in Eq. (3.15). Therefore, the re- 



sulting potential, scalar field and coupling function corresponding to the MCG shows 



the same dependence as given in Eqs. (3.19 3.22), and it can be shown that all the 



equations of the MCG, at B — become the corresponding equations for the GCG. 
Additionally to the GCG, the MCG model reproduces the radiation dominated phase 



of the universe at high redshift, as can be seen from (3.24) for B = 1/3 and a « 1 
(neglecting the constant term). In general the MCG reproduces the barotopic fluid 
with constant equation of state (p ~ (7 1 /( 1 + a ) a - 3 ( 1 + B )) a t earlier epochs when the 
constant term may be neglected. Resumming, the non-minimally coupled scalar field 
also reproduces the dynamics of the modified Chaplygin gas cosmology, with well 
behaved scalar field, potential and coupling function. 

The time variation of the gravitational coupling 

We can meet the constraints on the current value and the time variation of the 
gravitational coupling [51] . by appropriately defining or constraining the constants C 



and 9q. The effective gravitational coupling from (2.5) is given by 

G 



Geff ~ 1 - (3 ' 25) 
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where we used k 2 = 8ttG and <fi 2 = H 2 <p'2 = H 2 9. In terms of x, the time variation 
of the gravitational coupling can be written as 

G pff 3£k 2 d 



G, 



T eff 
eff 



3£k 2 FH 2 9 dx 



(FH 2 6)H. 



(3.26) 



Replacing the product FH 2 9 from Eq. (3.2), and evaluating at the present time 
(x = 0), the Eq. (3.26) can be written as 

G e ff 3/(C, ft 
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(3.27) 
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where f(C,6 ) = ^K 2 d(FH 2 9)/dx, g(C,6) = £k 3 FH 2 6 valuated at x = and the 
parameters of the model appearing in H 2 have been fixed, so that the resulting ex- 



pression (3.27) depends on the constants of integration C and 0q. We can meet the 
restrictions imposed by the current observations on the value and the time varia- 
tion of the gravitational coupling [51], if f(C,8) and g(C,6) satisfy the constraints: 
f(C, 0) ~ and g(C,9) ~ (these restrictions are actually softer: f(C, 9 ) < 1CT 1 



and g(C,0o) < 10~ 5 ). Thus, for the Chaplygin gas solution given by (3.6), (3.10 3.13) 
with A = 0.7 (B = 0.3), under above restrictions we found C ~ 11 and k 2 6o ~ 0.9, 
which give the potential plotted in Fig. 1. Following the same procedure for the 
generalized Chaplygin gas, we found C ~ 25.7, k 2 8 q w 1.5 (for A = 0.5 (B = 0.5), 
a = 0.6) and C ~ 3.1, k 2 9q w 1.5 for the same values of (A, B) and a = 3, with the 
respective potentials plotted in fig. 2. 



4 Discussion 

The cosmological implications of the Chaplygin gas model have been intensively in- 
vestigated in recent literature. We considered the model of scalar field with kinetic 
terms coupled non-minimally to the scalar field and to the curvature, to give a dynam- 
ical description of the Chaplygin gas model of dark energy. We have found analytical 
expressions for the reconstructed scalar field and potentials that describe the stan- 
dard, the generalized and modified Chaplygin gas models of dark energy dark mater 
unification. Thanks to the presence of the coupling function F(0), we could impose a 
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restriction on the potential through Eq. (3.4), which allowed us to find the solutions 
that lead to dynamical description of the Chaplygin gas cosmology. The results show 
that the obtained potentials V(4>) decrease with the evolving scalar field <fi. From 



Eqs. ( 3.10|3.12 ) it follows that is an increasing function of the redshift and V is 
a decreasing function of the redshift, which means that the scalar potential is a de- 
creasing function of the scalar field. In Fig 1. we show the V(<f>) dependence for the 
standard Chaplygin gas in the redshift interval [0, 2}. The same behavior follows from 
Eqs. ( |3.19|3.21[ ) for the GCG as shown in Fig. 2. The runaway behavior of all the 
potentials is a desired property of dark energy potentials as the relevance of the scalar 
potential increases at the present epoch (low redshift), as showing in Figs. 1 and 2. 
Although the problem of perturbations has not been considered here, it should be 
noted that, despite of the equivalence found between the specific reconstructed scalar 
model and the Chaplygin gas, under perturbations of the scalar field (which induce 
additional perturbations of the metric through the kinetic coupling) it may be ex- 
pected that this equivalence could not go beyond perturbative corrections. 
In order to satisfy the current restrictions on the Newtonian coupling, we can use 
the freedom in the integrations constants C and 9q, in order to control the actual 
value of G and it's time variation. These conditions can be satisfied by imposing the 
inequalities f{C, 0„) < KT 1 and g(C, 0„) < 1(T 5 . 

The above results show that the scalar field model with derivative couplings to cur- 
vature considered here, provide a dynamical scenario to describe the Chaplygin gas 
cosmology. The wide variety of phenomenologically acceptable solutions [3"4"l I33| l3"o] , 
support the capability of this model to explain the current status of the accelerated 
expansion of the universe, through different cosmological scenarios. 
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